Large negative velocity gradients in Burgers turbulence 
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We consider ID Burgers equation driven by large-scale white-in-time random force. The tails 
of the velocity gradients probability distribution function (PDF) are analyzed by saddle-point ap- 
proximation in the path integral describing the velocity statistics. The structure of the saddle-point 
(instanton), that is velocity field configuration realizing the maximum of probability, is studied nu- 
merically in details. The numerical results allow us to find analytical solution for the long-time part 
of the instanton. Its careful analysis confirms the result of [Phys. Rev. Lett. 78, 1452 (1997)] based 
on short-time estimations that the left tail of PDF has the form InP(u^) oc — [u^ 
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PACS numbers: 47.27.Gs, 05.10.-a 
I. INTRODUCTION 

We consider the random forced Burgers equation 



(1) 



that describes weak ID acoustic perturbations in the ref- 
erence frame moving with the sound velocity |^ . The ex- 
ternal force (f> in this frame is generally short-correlated 
in time, so let us assume that 



{(f>ixi,ti)(l){x2,t2)) = 6{ti - t2)xixi - X2). 



(2) 



Then the statistics of can be thought Gaussian and 
therefore is completely characterized by (^. We are in- 
terested in turbulence with a large value of Reynolds 
number Re = (x(0)L'*)^/"^/i^, where L is the characteris- 
tic scale of the stirring force correlator x(x). This prob- 
lem was intensively studied during the last years ||^-|^ . 

The main feature of Burgers turbulence is the forma- 
tion of shock waves with large negative velocity gradient 
inside and small viscous width of the front. The posi- 
tive velocity gradients are decreased by the dynamics of 
Burgers equation due to self-advection of velocity. On 
the contrary the increasing of negative gradients could 
be stopped only by viscosity. The motion of shock waves 
leads to a strong intermittency, the PDF of velocity gra- 
dients V{ux) is strongly non-Gaussian. The one way to 
describe the intermittency is to study rare events with 
large fluctuations of velocity, that give the main contri- 
bution to the high momenta (w") or to the PDF tails. 

The right tail (positive large u^) of PDF In 7^(1*2,) oc 
—u^ was first found by Feigel'man ||^ for the problem of 
charge density wave in an impurity potential. Later it 
was recovered using operatorproduct expansion |^ (see 
also 1^), instanton calculus j^, minimizers approach 
and mapping closure 0. The left tail in the inviscid 
limit seems to be algebraic, probably V{ux) cx \ux\~'^^'^ 
1^ (see also [0^). Due to viscosity the very far left 
tail is stretched exponential: h\'P{ux) (x —i''^\ux/i'\^- 
The large negative gradients exist practically only in- 
side the shock waves. The maximal value of gradient 



is proportional to the square of the velocity jump on 
the shock wave: |ua;|max = (Am)^/8j'. Then roughly 
the tail of the shock wave amplitude PDF has the form 
In'Pshock(Au) cx -i/3~2'3(Au)2'3. The analysis of the in- 
stanton structure predicts the value (3 — 3/2 This 
prediction is consistent with the assumption, that the 
tails of T'shocklAw) should not depend on the viscousity 



We are interested in the statistics of large values of 
gradients ^ u^rns/L ^ (x(0)/L^)"^/^. The velocity 
field configurations u{x,t) that make a contribution to 
the probability P(a) of the equality Ma;(0, 0) = o have the 
gradient greater or equal to a somewhere. The probabil- 
ity T'{a) decays very fast while a increases, i.e. the contri- 
bution of events with gradient greater than a somewhere 
is highly suppressed. Then one believes that only some 
specific field configurations u{x,t,a) ("optimal fluctua- 
tions" ||l^ or instantons) make contribution to P(a) at 
large a > {u-c-axs/ L) Re. Under this assumption to calcu- 
late V{a) one should find this optimal field configuration 
u{x, t, a) and estimate the probability of its realization. 

All instantons of this type are posed at the far tail of 
the statistical weight of averaging ^[(j>{x,t)]. Indeed, to 
produce large fluctuation of u the stirring force (jj also 
should be large, and the probability of such fiuctuation 
is low. The weight may not contain a large param- 
eter, but it should have fast tails, e.g. exponential ones. 
Then the concurrence between statistical weight and the 
value of calculated quantity makes the contributing real- 
izations of 4>{x, t) rather determined. This approach was 
introduced by Lifshitz . Later it was applyed to deter- 
mine high order correlation functions in field theory [Tl| ] 
and in the systems of hydrodynamic type: simultaneous 
(see, e.g. ]l2|,|jl3)) and non-simultaneous |l^] ones. 

The paper is organized as follows. In Sec. || we de- 
rive the equations for the instanton. Sec. |l| is devoted 
to the detailed description of our scheme of numerical 



calculations. In Sec. IV we discuss the numerical results 
and describe the behavior of the solution of instanton 
equations at large times. 
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II. SADDLE-POINT APPROXIMATION 

The velocity gradients PDF V{a) can be written as the 
path integral 

ioo 

= J VuVp J dJ^ exp(^- S + 4i/2jF(u:r(0, 0) - a)) , (3) 



where the effective action S has the form [y_5 12 



5=i J dt J dxidx2p{xi,t)x{xi ~ X2)p{x2,t) 

— oo 



- i dt dxp(ut + uux - vuxx)- (4) 



The integration over T gives rise to 5(u2;(0,0) — a), and 
the factor AiP was chosen for our convenience. Note that 
if retarded regularization of the path integral is used 
then J VuDp exp(— 5) = 1 and we have no normalizing 
M-dependent denominators in (|^). One can find some 
analogies between the appearance of the second field p 
and technique that was developed by Keldysh [|l6| for 
nonequilibrium dynamics description. 

We are interested in the tails of PDF 7^ (a), i.e. the 
parameter a in the integral (H) is large. The asymptotics 
ofV{a) at large \a\ > (x(0)L)2/3/ u is determined by the 
saddle-point configuration of fields u{x,t), p(x,t) (and 
also parameter T), near which the variation of the inte- 
grand is equal to zero [|l2|. The saddle-point configura- 
tion (sometimes called classical trajectory or instanton) 
is governed by the following equations 

ut + uUx " vuxx = -ix* P, (5) 
Pt+upx + vpxx^'^iv^J^5{t)6'{x). (6) 

where % * p is the convolution 

(X * P){x) = j dx' x{x - x')p{x'). (7) 

The solution should satisfy boundary conditions 
lim u{x,t) = 0, lim p{x,t) = 0, 
lim u{x,t) = 0, lim p{x,t)=0. (8) 

\x\ — *oo \x\ — »oo 

The value of !F is tuned in such a way that the condition 
Ux(0,0) = a holds. The quantity is a Lagrange mul- 
tiplier for finding the extremum of S with the condition 
UxiO,0) = a. 

The equation for p should be solved moving back in 
time because of the signs at pt and Pxx in the instanton 



equation (^). The convolution —ix*p is the optimal con- 
figuration of external force </> that produces large negative 
gradient. 

In what follows we will measure the length in L units, 
i.e. we set L = 1. Rescaling the time t and fields u, 
p one can exclude the parameter v from the instanton 
equations: 

t = T/2iy, u = 2i^U, p = Aiiy'^P, a 2i^A, (9) 

Ut + UUx ~ ^^Uxx ^ j dx' x{x~x')P{x'), (10) 

Pt + UPx + \Pxx^TS{T)5'{x), (11) 

at T = one has [/j;(0, 0) = A. The only parameter 
in the instanton equations is ^ = a/2v. Note that the 
steady-state kink solution of Burgers equation with the 
negative gradient a is 



= — ■\/2iy|a| tanh (^\J\a\/2vx 



(12) 



Thus the physical meaning of |A| is the square of the ra- 
tio of pumping scale L — 1 and the kink width Wkink = 

The effective action iSoxtr at the instanton that gives 
the right exponent: 



InP(a) ~ -5cxtr(a), 



(13) 



is equal to 



5extr = / / dxidX2p{xi,t)x{xi - X2)p{x2,t) 



u 

4zy3y" dT j dxidx2 P{xi,T)x{xi - X2)P{x2,T). (14) 



The freedom of rescaling the fields u, p and the time t 
with appropriate change of v gives us the following rela- 
tion: 



5extr(a) = 8z/35(a/2zy) = (2i/)35(A), 



(15) 



with the function S{A) to be determined. One can prove 
by straitforward calculation the following relation be- 
tween functions T{A) and S{A): 



T{A) 



dS{A) 
dA ' 



(16) 



The relations of such sort are well-known in classical 
mechanics; here A and J- are conjugate variables, and 
saddle-point configuration is the trajectory of extremal 
action. 

The instanton equations ( p^ , pl] ) are Hamiltonian: 
Ut{x,T)^~-^^, Pt{x,T)^-JP—, (17) 



5P{x,Ty ' ' 5U{x,T) 

j dxP{uUx-\Uxx-\x*P)- (18) 
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The Hamiltonian 7i is the integral of motion, i.e. 
dJi/dT — 0. Since both U and P tend to zero at 
T — !• — cx) we have 7i = 0. From the instanton equa- 
tions and the condition 7i = we get 



S=^ + \ j dTdxP^U'^^ 



i / dTdxP.^U^. 



(19) 



The last term is due to viscousity. At the right tail 
it is unimportant, and we have dS/dA = 3S'/A, i.e. 
S . At the viscous left tail its contribution to the 
action is of the same order as other terms. If In'P(a) 
is a powerlike function: ln7-'(a) oc \a\^ , then one has 
/ dTdxP^U^ =2(3-/3)5'. 

The high momenta can be calculated by the instanton 
method in a following way. Because a^V{a) is a narrow 
function for large n, and only narrow velocity interval, 
which position depends on n, contributes to (a"). The 
position of this interval is exactly the saddle-point in the 
integral (a") oc / da a" exp(— 5oxtr(a)) (see (|l^)), that 
satisfies the equation 



d5cxtr(a) 
da 



dS(A) 
dA ' 



(20) 



Combining it with T = n/^v'^A we again get ([l^). To get 
the instanton equations for the average (a") one should 
only substitute T in ( [lT| ) for n/div^A. Then the instanton 
equations become the same as in Q. 

One also should consider fluctuations near the instan- 
ton as a background. The way how the fluctuations can 
be taken into account is unknown yet but their influence 
to \nV{a) due to their phase volume is small in com- 
parison with iSoxtr while a 3> (urms/^) ri-e. At smaller 
gradients the fluctuations essentially change the answer 
and we have the algebraic tail . 



III. NUMERICAL CALCULATIONS 

The preliminary calculations that were made in x, T 
variables have shown that the width of the instanton 
equations solution grows with |r| and is proportional 
to |T|^/^, while its amplitude is proportional to |T|~^/^. 
To avoid the necessity of treating simultaneously narrow 
structure at small T and wide one at large T we used the 
following variables: 



- T, T = To (1 



U = U/^To - T, P^P/^To-T, 



(21) 
(22) 



where Tq is some constant of the order of unity. The 
instanton equations in these variables take the form 



Ur + ^{^U- U^) + UU^ = x(t) * P, 



(23) 



Pr + ki^P + Pi). + UP^ = T5{t)6'{OIV%. (24) 



where x(C, t) = (To — T)'^/'^x{^)- The boundary condi- 
tions for U, P are analogous to (||). 

Let us describe now the general structure of the nu- 
merical scheme that finds the solution of our boundary 
problem. 

The diffusion terms U^^, in instanton equations 
(^,^J) have opposite signs. If one considers these equa- 
tions as two linked Cauchy problems, then the natural 
direction of time in (^3|) is positive, while in ( p4[ ) the 
direction is negative. Assume that at a given value of 
T the approximate solution t/oid(Cj''") is known. Let us 
try to make it closer to the true solution of the prob- 
lem. For this purpose let us solve the Cauchy problem 
for ( |2^ ) starting from r — -|-0 and moving up to large 
enough Tmin < 0. Then using P that we have got in the 
previous step we solve the Cauchy problem for ( p3| ) mov- 
ing from r = Tmin up to r = 0. As a result we get the 
new values C/now(C; '^)- Further we will use the sign / for 
the mapping Uo\d U-acw The stationary point U of 
the mapping / and the corresponding function P are the 
desired solution of ( p3| , p4| ). 

The numerical experiments have shown that iterations 



U 



(25) 



converge ii T > T^, ~ —0.96. While T < the simple 
iterations ( p5| ) are divergent. 

Curve 1 at the Fig. ^a) shows how the value of the 
gradient A = 9C/(^, r)/9^|^=o,T=o depends on the num- 
ber of iteration for = — 2. It can be seen that the sta- 
tionary point of / is unstable, but the mapping /[/[C^]] 
has two stable stationary points. The stability properties 
of the iterations are determined by the spectrum of the 
linearization K of / near a stationary point: 



/[c7 + y] = ij ^kv 



(26) 



The iteration process is convergent if the modulus of all 
eigenvalues of the linear part K is less than 1. The period 
doubling indicates that while J- passes through one 
of i^'s eigenvalues passes through —1 [0. Let us denote 
this eigenvalue as A. 

This knowledge allows us to construct the new map- 
ping with stable stationary point that coincides with one 
of the mapping /. Let us proceed the following iterations: 



[7(^+1) = c/[J7(*'] + (1 - c);7''\ 



(27) 



where c is some constant. It is easy to check that the sta- 
tionary points of the mappings / and /c do coincide. The 
linear part of /c is equal to K^. = cK+ (1 — c)l, its eigen- 
value corresponding to unstable A is Ac = cA -I- (1 — c). 
If we take the value of c inside the interval < c < 
2/(1 - A) < 1, then |Ac| < 1, and the iterations (|2^ are 
convergent. 

Since we don't know A a priori, the value of c that 
provides the convergence of iterations, was determined 
from experiment. Fig. ||(a) illustrates the influence of 
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the decreasing of c on the dependence of A vs. number 
of iteration. In the final version of the computer code the 
value of c was changed in an adaptive way: each time the 
value of 1^1 was decreased after an iteration, c was mul- 
tiplied by 0.9. One can compare from the Fig. |l|(b) the 
iterations run for c = 0.1, c = 0.05 and for adaptive de- 
creasing of c, all three for = — 2. The initial and final 
values of c in the case of adaptive change were equal to 
0.1 and 0.1 • 0.9" = 0.05905..., respectively. It is clear 
that for the two last cases the iterations converge to the 
same solution. 




10 20 30 40 50 
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-100 
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that in principle are posed at the infinity were consid- 
ered as they are posed at (large enough) ^,nax and T^-m- 
Typical values used: ^max = 10 and Taun ~ —30. The 
grid had uniform mesh intervals in ^, typical number of 
grid sites along ^ axe was eqial to 1024. 

The first calculations had shown that the solution 
changes rapidly in the vicinity of r = 0, while at large 
T it varies slowly. Because of this we used nonuniform 
grid for variable r. The time step was smaller inside 
the interval ti < r < 0, typical value used: t\ = —0.4. 
The number of grid sites inside this interval varied from 
2000 to 4000. The computer power we had limited the 
number of all time steps in the grid by 5000. During all 
calculations we used Tg = 1. 



B. Cauchy problem for P 

The equation ( p^ should be solved backward in time. 
The source term in right hand side of the equation (|2^ ) 
provides us the initial condition P{t — —0) (x 5'{S). This 
means that at small times the field P is localized in a very 
narrow interval centered at ^ = 0. Such initial condition 
can not be accurately discretized, so we used another way 
to represent P at small times. 

For small r the field P is very narrow. At its support 
we can approximate the velocity J7 by a linear profile: 
U — A{t)^. The evolution of P in such a velocity field is 
described by the derivative of Gaussian contour: 



/'amp(-O) = 



e^amp(r) 

^27rTo7^3(T; 
1, D{-Q). 



: exp 



0. 



e 



2D{t) 




1000 2000 3000 

Number of iteration 



4000 



5000 



D{t) = / dr' exp / dr" {2A{t") 



dr'(2i(r') + 1/2) 



-Pamp(r) = exp 



(28) 



(29) 



(30) 



FIG. 1. The gradient A as a function of number of iter- 
ation. At Fig. hia) the dots corresponding to one run are 
joined by a line for convenience. The curves are: (a) simple 
iterations according to ( p^ ) (curve 1) and according to ( |2^ ) 
with c = 0.1 (curve 2); (b) c = 0.1 (curve 1, its thickness is de- 
termined by the amplitude of serrated oscillations), c = 0.05 
(curve 2), and adaptive decreasing of c during the run, the 
initial value is c = 0.1, c was multiplied by 0.9 four times 
(curve i); T = —2. 



A. Grid parameters 

The solution of Cauchy problem is found numerically 
using the method of finite differencies. The grid covers 
the rectangular domain < ^ < ^max, Tmin < r < 0. In 
numerical calculations some of the boundary conditions 



We use such a representation for P for ro < t < 0. The 
value of To is chosen in such a way, that the velocity field 
U is still linear at the width of P. From the other hand, P 
at T = To already becomes wide in comparison with mesh 
interval A^. Typical value of tq used: tq = —1/1500. For 
times r < tq the solution is found by fully implicit scheme 



P, 



n+l 



pn 



At 



2 ' 



D 



2P" 



pn+l 



pn+l 

-1 



Ae 

pn+l _ pn + l 
n+lfi ^i-l 



0, (31) 



where ^ , I?" are equal to 



= 



0.5(r(e„T„)±|r(C^,r„ 
1 

l + 0.5|r(C„r„)|A^- 



(32) 
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The function r(^, t) is expressed via velocity field: 
r(e,T) =0.5^ + ?7(e,T). 

Here > 0, At < — mesh intervals, and 
Tn — site coordinates. The numerical scheme used is 
monotonous and stable, it is of the first order of accu- 
racy in At and of the second order in A^ . 



IV. VISCOUS INSTANTON 

In this section we represent the results of our calcu- 
lations, that show the structure of the instanton and its 
change with \J-\. The minimal value of J- at which the 
reliable results in numerics were obtained is = —2. 



C. Cauchy problem for If 

At this stage we use the initial condition U{t — Tmin) = 
0. The viscousity, source and self-advection terms are 
treated by splitting technique At each time step it 
is first calculated the change of U due to the source, then 
— to the viscousity, and at the last — to the nonlinearity. 

From calculated already grid layer ?7" — U{t — Tn) 
the next layer U"^-^ is found in a following order: first, 
the equation Ur — %(''') * P: is solved: 



U, 



ri+l 



ix{r) * P)n. 



(33) 



The convolution is calculated as the result of inverse 
fast Fourier transform (FFT) acting on the product of x's 
and P's FFT images. The external force correlator dur- 
ing all calculations was equal to x(a;) = (1 — x'^)e~^ = 
—(Pe~^ 1"^ jdx^ . The numbers C/f'''^ that are found in 
such a way are not a final solution for the layer [Z"**"^, 
since only the source term has been taken into account 
yet. We use them as an input at the next step, we will 
denote them as J7" (note that they do not coincide with 
C/fin||)). 

Next, the viscousity and linear part of advection are 
taken into account, according to the equation Vt ~\- 
\{iJJ — U^^ = 0. The fully imphcit scheme was used 
analogously to (pl^) : 



ur 



At 



2 ' 



D, 



U. 



71+1 

i+1 



+ 



n+l 



n+l 
1 



AC 



(34) 



here A = 1/(1 + 0.2^^iA^). Again, the numbers C/"+^ 
do not form a final solution, and we send them to the 
next step with a C/" notation. 

The nonlinear part of the equation Ut + UU^ = was 
solved by explicit conservative scheme pO|: 



At 



2A£, 



= 0, 



(35) 



that finally gives us the next layer of the velocity field 
Ijn+i^ This scheme is of the first order of accuracy in At 
and of the second one in A^. 




FIG. 2. The level curves if U{^,t) for = -0.9 (a), 
= —1.1 (b), and T = —2 (c). The values of levels can be 

calculated from the two given levels according to arithmetic 

progression law. 

The general features of the instanton structure change 
with can be obtained from Fig. |^ that shows the level 
curves of U{£,,t) for three values of J^. Since U and 
P are the odd functions of C, we draw only the region 
where ^ > 0. The calculations were done in a rectangu- 
lar < ^ < 10, —30 < T < 0, whose dimensions are a bit 
larger than it is shown at the figure. 
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One can see that the instanton hfe-time and the max- 
imal value of \U\ rapidly increase with the growth of 
\J-\. The growth leads to the deformation of level curves 
near r = because of the influence of nonlinearity, that 
is weak at = —0.9 (Fig. |^(a)) and very strong at 
(Fig.ic)). 

A. Structure 

Detailed analysis of the instanton solution based on the 
results of numerical calculations allows us to distinguish 
five different regimes in the instanton time evolution. Be- 
low we discuss them consequently from t — to t — -co. 




-2 




FIG. 3. The kink in the velocity field U{^, r) at r = 0; 



considering evoWtton back in time they are expanded by 
it since Ux\x=q < 0. The time needed for the expansion 
is equal to - L/U - l/v^PT- 

During the 3rd-5th regimes the width of U and P field 
is much greater than L = 1. Then it is naturally to sub- 
stitute x{^) by —X2S"{x). The instanton equations take 
the following simple form: 

Ut + UUx + ^Vxx=0, Vt + UVx + ^Uxx = 0, (36) 

where V — 2x2P — U, X2 = ^5/ dxx^xix)- While 
moving to large negative time numerical solution has a 
tendency to fall in U = V (see Fig. ^, curves r = —4 
and T = —8). Then the equations for U, V reduce to the 
Burgers equations with the evolution back in time. 




FIG. 4. Comparison of C/(^,t) (solid curves) and X'iPi^y''') 
for three values of r; T = —2. 



The first regime consist in the viscous smearing of 
the field p up to the scale of the kink width Wkink — 
^/2vl\a\ = \I^J\A\ (see (H), Fig. |). Since the viscousity 
plays the crucial role in this regime, we will also use di- 
mensional variables. kX t = —0 we have p{x,t) cx S'{x), 
and the width of the kink in the velocity profile is equal 
to Wkink = •\/2i//|a| (see (p^)). Since p is very narrow, 
the viscousity dominates the evolution. The width of 
p obeys the diffusion law and equals to y/2iy\t\. These 
two widths become comparable at time t = —l/\a\, or 
T = —1/\A\. This means that during the whole time of 
smearing of p by viscousity the width of the kink is of the 
order of u>kink- Indeed, if the shape of the velocity profile 
deviates from steady-state kink solution (|lj), then the 
change of the kink width during this time period would 
be of the order of ut ~ •\/2z^|a|/|a| = y^2z7/|a| = Wkink. 
At this regime the source term x * -P in instanton equa- 
tion ( [lO| ) is unimportant. When the width of p becomes 
of the order of Wkink the rate of expansion of p due to the 
velocity gradient becomes comparable with the rate due 
to viscousity (such a balance determines the width of the 
kink). 

The next (second) regime was exhaustively studied in 
It consists in dilation of fields P up to the pump 
scale L = 1. The fields are advected by velocity C/, and 



Such a substitution is not always possible. During such 
an evolution the shock waves occure (see, e.g., the level 
curves at Fig. |^(c) near t = —5, Fig. ^, curves 1 and 2). 
For transition to equations (^) to be valid the width of 
these shock waves should be larger than L = 1. Other- 
wise the substitution of x by 5" is not valid. The width of 
the shock wave in dimensionless variables is greater than 
1 only if its height is smaller than 1. However, right after 
the 2nd regime when the width of U and P fields becomes 
greater than the pumping force correlation length L = 1 
the amplitude of the velocity field U is of the order of 
yJ\A\ ^ 1. The amplitude of U becames small only at 
very large times (it is shown below that such a crossover 
happens at T ~ — ^|yl|). It means that there is an in- 
termediate regime that goes after the 2nd one, where the 
substitution of x by 5" is inapplicable. In this (third) 
regime the fields U and P arc smooth functions in the 
interval wider than 1. At the ends of this interval they 
contain shock waves — the value of U and P rapidly goes 
to zero (as it is shown schematically in Fig. ||). We will 
use the word "shock" for these structures at the ends of 
the interval, while for narrow structure in the velocity 
field U near a: = 0, T = we will use the word "kink" . 

Now we will consider the structure of the shocks in de- 
tail. Let us denote the height of the shocks in U{x, T) 
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and P(x,T) fields by Hu{T) and Hp{T), respectively. 
Their position we will denote as Xshock(r) (xshock > 0, 
shocks are posed aX x — ixghock)- 




FIG. 5. Schematic representation of U and P for the 3rd 
regime, T — > — oo. 




FIG. 6. Shock structure in the 3rd regime, T — oo. 

In the third regime the width of the shocks in the 
field P is determined by the competition of squeezing 
them by the velocity U and spreading by viscousity (see 
Fig. ^). Since we are in a strongly nonlinear regime 
the viscousity is weak, and the width of P's shocks is 
very small. The shocks in P are stationed at the cen- 
ter of U's shocks, i.e. in almost linear velocity profile 
with the gradient Ux ^ Hu- The width of P's shocks 
can be estimated as l/^/Hu ^ L. Then the good ap- 
proximation for P{x,T) near the shock x ~ a^shock is 
P{x,T) ~ ~Hp{T)6{xshockiT) - x), where e{x) is step- 
function. During the evolution forward in time the shocks 
in U do not break down because of source term x * P. 
The source prevents the destructive effect of advection 
term UUx- The shocks in P should carry the [/'s shocks 
of the height Hjj, i.e. of strength UUx ~ H^. Thus we 
should have Hp cx H^j in this regime. Now we show it 
more carefully. 



Let us write the instanton equation (|0|) in the refer- 
ence frame of the shock near the point x = Xshock (see 
Fig. |6|). We have two contributions to time derivative Ut- 
from the growth of Hu in time (of order Hjj /T) and from 
the motion of the shock (of order H^). Neglecting the 
first one we write the following equation for the velocity 
f/(x,T): 

1 2 

-{U{x) - f7(xshock))^ = -HpX'{x - Xshock), (37) 

where the new function X{x) is determined by the equa- 
tion x{x) = '^^"{x) with the condition X ^ Q with 
X — > ±oo. Integrating this equation once we obtain 



{U{x) - ?7(a;shock))' 



2Hp{X(0)-X{x~x,Y.ock)). 

(38) 



Since C/(a;shock) = -Hu/2 we have Hp = i/2/8X(0). 

The next step consists in finding the solution of (10 1^) 
between the shocks posed aX x = ixghock considering the 
fields U , P as smooth ones and using the boundary con- 
dition 

P(±a;,hock T 0, T) = tU^{x,^,,cV - 0, T)/8X(0). (39) 

The instanton equations (|l^,|ri|) take the form 

Ut + UUx = 0, Pt + UPx = 0. (40) 

Here we approximate the shocks as jump discontinuities, 
and the condition (|3^ ) relates the heights of the jumps 
(see Fig. H). Here the diffusion terms and the term x * P 
are omitted. One can check that they are negligible since 
the characteristic x-scale of the solution is large enough. 

The equations (^) can be integrated by characteristics 
(or Lagrangian trajectories). The velocity of the shocks 
is equal to ztHu/2, i.e. all the trajectories disappear at 
the shocks (if we consider the evolution back in time). 
The value of U (or P) is conserved in time if we follow 
the Lagrangian trajectory. This means that the relation 
P = U\U\/8X{0) holds everywhere between the shocks. 

Due to self-advection the velocity field U becames 
more and more linear as a function of x while |r| in- 
creases. This happens at the border between the 2nd 
and the 3rd regimes. In the third regime we can take 
U{x,T) = x/T between the shocks. The field P is equal 
to P{x,T) = -x\x\/8X{0)T'^. The velocity U simply 
squeezes or expands the field P without changing its 
shape. Since Hp ^ H^, the field P should have the 
same scaling as J7'^, i.e. P (x x^. The concave form of 
P (as one can see ^ = —2 is not yet good enough for a 
clear picture) is shown in Fig. ^. 

Let us determine the time dependence of Xshock- We 
have 



^^shock 1 ^rTi\ a^shock 

-^^--Hu{T)-^^^. 



(41) 



Solving this equation we get Xshock(P) = B\/—T. Since 
xshock - 1 at |r| - l/^/\A\, we get B - \A\^^^. The 
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shocks heights are equal to Hu{T) = B/^f^, Hp(T) = 
BV8X(0)|T|. The width of the shocks in P field is of 
the order of 1/ y/Hu- The shocks in velocity field U have 
the width of the order of 1, since U is pumped by x * P. 




U{^) (solid curve) and x^PiO (dashed curve) pro- 
= -0.4. :f= -2. 



FIG. 7. 
files at r 



At large time T ^ — y'I^I the height Hjj (and con- 
sequently the width of the shocks in P) becames of the 
order of 1. This indicates the end of the third regime and 
the beginning of the fourth. Going further back in time 
we finally enter the domain of validity of the equations 
(^). The solution falls into U = V. Again, the solution 
has two shocks between which it is a smooth function. 
The shock position satisfies Xshock oc ^/—T, between the 
shocks U = V = x/T holds. This regime exactly corre- 
sponds to self-similar solution u{x, t) = 9{t — Cx^) x/t of 
inviscid Burgers equation ut -f uux — • u^x = (here C 
is a parameter). 

At the far tail of the instanton due to viscous dissipa- 
tion the solution U = V transformes to the derivative of 
Gaussian contour — the self-similar solution of the dif- 
fusion equation (see Fig. ^, r = —16). During the fifth 
regime the advection term UUx becomes irrelevant. In ^, 
T variables the solution tends to U (x ^exp(r/2 — C^/2), 
that was observed in numerical calculations. 

During the 4th and the 5th regimes the amplitude of 
the velocity field U is less than unity. It means that the 
amplitude of U is lower than the order of the typical sta- 
tistical fluctuations, and the saddle-point approximation 
is meaningless there. The typical events that demon- 
strate large negative gradients start from some velocity 
configuration U(x) with the amplitude of the order of 
unity and are governed by the 3rd regime first. The ac- 
tion on these events and their further evolution almost 
do not depend on the initial velocity field U{x), and the 
dependence oiV{a) on a remains unaltered by the averag- 
ing over all possible configurations U (x) . We considered 
the 4th and the 5th regimes since they are the parts of 
the whole solution of our nonlinear boundary problem. 

Let us now run the whole evolution forward in time. 
At the beginning (5th regime) the field U is pumped by 



very wide P. The pumping force is proportional to P^x- 
During the 4th regime the source is localised at shocks 
in P that leads to a formation of shocks in U. The J7's 
shocks want to break down because of self-advection, but 
the source term x* P keeps them going. When the grow- 
ing height of P becomes larger than unity the shocks 
in P become narrow. The balance between the terms 
UUx and x * -P in (0) changes a little, that results in 
the change of the form of U's shocks — it is determined 
by the shape of xi^) now. The distance between the 
P's shocks decreases in time and eventually it becomes 
comparable with unity. After this P becomes even more 
narrow and the efficiency of the source term begins to 
fall down. The self-advection of the velocity destroyes 
the shocks and leads to a formation of the kink at 2: = 0, 
while P transformes to 6'{x). The kink shape at JF = —2 
is shown in Fig. ^ Schematically time evolution of the 
instanton is illustrated in Fig. |[ 



1st 



2nd 



3rd 



4th 5th 



U 
P 



-0.5 

T = n 



0.5 

T = Tb 



ln|T| 
In 1^1 



FIG. 8. Schematic representation of the instanton time 
structure. 



B. Action 

One can present the action S{A) in the form of S* = 
/"^ dTs{T) like in expression (|l|). For = -2 the ac- 
tion density s(T) that was obtained from numerical cal- 
culations is shown in Fig. ^. While T < ^ -l/^/pJ 
the convolution % * p is localized at shocks, so s{T) ~ 
H'I,{T) - B-^/T^. The maximum of s(T) is posed at 
T ^ T^,. Further increasing T leads to the decreas- 
ing of the density s{T) because P{x) becomes more and 
more narrow without an adequate growth of its ampli- 
tude. This region of small time T > was studied in 
Q. It was shown that the contribution St>t, to the ex- 
tremal action from this interval is of order and the 
main contribution to the action from it comes from the 
region T ^ T<, — the border between the 2nd and the 
3rd regimes. Exactly these two regimes determine the 
optimal configuration of noise providing the event with 
large negative gradient. 

The contribution of region of time — yJ\A\ = Tb < T < 
r» = — 1/-\/|^| (3rd regime) to the extremal action S{A) 
can be estimated as 



dTHl{T) 



-B^/T^ ^ \A\ 



3/2 



(42) 



Tb 
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Note that the value of St<t, is again cumulated from 
the region T ^ T^,. The crucial point is that the contri- 
bution to the action S{A) from the tail of the instanton 
(or large time T < T,) is finite, i.e. the integral (^ ) con- 
verges (the addition to the action from interval T < Tb 
is negligible). Also this contribution is not dominant, i.e. 
it is not much greater than the contribution of the or- 
der from small times (T > T*). It means that in 
our case the instanton is localized enough in time. Its 
long-time dynamics does not destroy the fact that it is 
the main fluctuation determining the statistics of large 
negative gradients. 




i\ s 
- 500 



-0.5 Th, Ol r 

FIG. 9. The action density s(r) as a function of "time" r. 
J^ = -2. 



tial condition for U during iterations. It turned out, that 
e.g. for = —2 we get the following values of TAjS and 
A with different value of t^u\- 



'^min 


TAjS 


A 


-4 


1.441 


-148.1 


-30 


1.437 


-102.6 



Although, when calculating with Tmin = —30, the con- 
dition C/(^, -4) = X2P{i, -4) holds within 15% (as one 
can see from Fig. we prefer to use the grid shorter in 
time (rniin = —4) to have smaller time step. The value 
of A strongly depends on time step. It was observed 
in numerical experiment. Such sensitivity is characteris- 
tic of Burgers equation also. Although the calculations 
with small |rinin| give us worse accuracy at the tail of the 
instanton, the smallness of the time step allows us accu- 
rately describe the main part of the instanton where the 
nonlinearity level is high. 

One can see the cubic asymptotics 5 cx at A > 0. 
The instanton structure for A > that was described in 
was confirmed by our numerical calculations. The 
case A < corresponding to the PDF's left tail is more 
complicated. The function TAjS has minimal value at 
A ~ —12. At further decrease of A it starts to grow and 
finally tends to the value 3/2. In this case the coefficient 



V. CONCLUSION 



TAIS 



° o oooof 

1 



o o °°° 



3/2 



-100 

FIG. 10. TA/S = d{lnS)/d{\nA) (see (^) as a function 
of gradient A. 

At the Fig. 0the function d{ln S) / d{ln A) = J^A/S 
that was obtained from numerical calculations is shown. 
We used different grid parameters during calculations 
for instanton structure and for this figure. Here we 
used Tmin — —4 with boundary condition C/(C,Tniin) = 
X2^(C; ''niin)- This boundary condition was used as ini- 



We have examined the remote left tail of the velocity 
gradients PDF V{ux) in Burgers forced turbulence. The 
possibility of direct numerical solving of instanton equa- 
tions by iterations is demonstrated. Numerical calcula- 
tions and the analysis of the instanton behavior at the 
time large compared with its lifetime ^ l/V ^\ux\ with 
the solution at small time from H] show that \a.V{ux) (x 

-h,,j3/2 
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